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In quantum held theories with topological sectors, a non-perturbative quan¬ 
tity of interest is the topological susceptibility Xt- In principle it seems 
straightforward to measure Xt by means of Monte Garlo simulations. How¬ 
ever, for local update algorithms and hne lattice spacings, this tends to be 
difficult, since the Monte Garlo history rarely changes the topological sec¬ 
tor. Here we test a method to measure Xt even if data from only one sector 
are available. It is based on the topological charges in sub-volumes, which 
we denote as slabs. Assuming a Gaussian distribution of these charges, this 
method enables the evaluation of Xt, as we demonstrate with numerical re¬ 
sults for non-linear a-models. 
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1 The topological susceptibility Xt 

There are a number of models in quantum field theory, which have the prop¬ 
erty that the conhgurations occur in distinct topological sectors; each sector 
is characterised by a topological charge Q G Z. This refers either to inhnite 
volume and conhgurations with hnite actions, or to hnite volume with pe¬ 
riodic boundary conditions. Here we consider the latter setting, with some 
volume V in Euclidean space. 

The models with this property include in particular 4d SU{N) Yang- 
Mills theories for all Y > 2. Fermions may be present as well, so this class of 
models encompasses QCD. In that case, the quenched value of the topological 
susceptibility Xt has a prominent application in the explanation of the heavy 
mass of the pseudo-scalar 77 '-meson [T]. 

The measurement of Xt is a non-perturbative issue, hence lattice simu¬ 
lations are the appropriate method. If the Monte Carlo history changes Q 
frequently, this measurement is straightforward. However, for most of the 
popular algorithms, including local update algorithms, as well as the Hybrid 
Monte Carlo algorithm (which is standard in QCD with dynamical quarks), 
the auto-correlation time with respect to Q tends to be very long, in partic¬ 
ular on hne lattices^ This problem is getting even worse in the presence of 
chiral fermions. 

Recently, several indirect methods to measure Xt nevertheless have been 
suggested and tested mu- Here we address a different approach for this 
purpose, which was hrst sketched in Ref. [^, but which has not been ex¬ 
plored ever since (though a similar approach was studied last year [ 8 ]). It is 
described in Section 2 . Sections 3 and 4 give results for the Id 0(2) and the 
2d 0(3) model, respectively. Our conclusions are discussed in Section 5, and 
an appendix is devoted to analytical considerations about artifacts in Xt- 

^In general, the topological sectors are well-defined only in the continuum limit (there 
are exceptions for lattice actions with a constraint that only admits very smooth con¬ 
figurations US]). So in general transitions are only enabled by lattice artifacts, and by 
discrete jumps of the algorithm. 
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2 Evaluating Xt from slabs 


We assume parity symmetry to hold, which implies (Q) = 0. We further 
assume the topological charges to obey a Gaussian probability distribution!^ 


p{Q) oc exp {-Q^/{2xtV)) , Xt = ^ (Q^) ■ (2.1) 


The idea of the method that we are going to explore, is to divide the 
periodic volume V into sub-volumes, which we denote as slabs, and to extract 
Xt from the fluctuations of the topological charge within these slabs. This 
has the potential of providing the result even based on configurations of a 
single topological sector (with respect to the entire volume V). 

We consider just two slabs, of sizes xV and {l—x)V, where 0 < x < l]^ At 
hxed charge Q, we denote the topological charge contribution in the hrst slab 
as g G IR; it is obtained by integrating the topological charge density over the 
slab volume. In general it is not an integer, because not all slab boundaries 
are periodic. Thus we obtain the probability for the charge contributions q 
and Q — q in the two slabs. 


vM) ■ V 2 {,Q - q)\x,v oc exp 

= exp 

oc exp 


1 r g' ^ 
2xtV U(1 - x) 


{Q-q? \ 

2XtG(l -x)/ 



2xtV x(l - x)7 


( 2 . 2 ) 


where we defined q' := q — xQ. 

Eq. fl2.2p implies (g) = xQ, and therefore 


(,'") = ( 9 =) 


(2.3) 


Thus, if we measure (g^) at a set of x values in a fixed sector with topological 
charge Q, we can fit the results for (g'^) to the parabola XtVx{l — x), which 
is predicted by eqs. fl2.2p and fl2.3p . 

^The assumption of a Gaussian distribution of the topological charges appears natural 
in light of instanton gas models. In case of the quantum rotor in infinite volume, it can 
be demonstrated analytically [S]. Simulations confirm that it also holds — at least to a 
good approximation — in the 2d 0(3) model [6] (cf. Section 4), and in QCD [TO] , 

^The extension of this method to a larger number of slabs is straightforward, but hardly 
promising, since the slab volumes should not be too small. 
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The simplest case is the topologically trivial sector, Q = 0, where (g^) (x) = 
{q'‘^){x) is given by a parabola, which takes its maximum at x = 1/2, with 
a value of (Q^)/4. In the topologically charged sectors, (g'^)(x) still has the 
same shape, whereas (g^)(x) is a parabola that connects (g^)(x = 0) = 0 
with {g^)(x = 1) = Q^. In any sector, the measured data for (g^)(x) can 
be htted to the expected parabola. The susceptibility Xt is the only htting 
parameter, which is evaluated in this manner. 


3 Results for the quantum rotor 


We hrst consider a quantum mechanical scalar particle moving freely on a 
unit circle. We deal with periodic boundary conditions in Euclidean time. 
Then there is an integer winding number for each trajectory, which charac¬ 
terises the topological sectors of this model. It is also denoted as the Id 0(2) 
model, or Id XY model, and it is related to 2d U(l) gauge theory [TT] . 

Now we assume Euclidean time to be uniformly discretised, and on each 
site t = 1,2... L there is an angular variable 4>t (with 4’t+L = We apply 
the geometric dehnition of the topological charge Q [12] . 

1 ^ 

Q[(P] = ^ 271 e (-TT, 7l] , (3.1) 

t=l 

where we dehne the modulo function such that Acpt is the shortest arc length 
connecting nearest neighbour angular variables. 

In our numerical studies we simulated three lattice actions: the standard 
action, the Manton action [13] and the constraint action 00 


‘S'standard [0] 

‘S'Manton [0] 
‘S'constraint [0] 


^ (l - cos(A0i 


t=l 

t=i 

0 < S \/t 

-l-cxo otherwise. 


(3.2) 


'^The constraint action is a special case of a topological lattice action mMm, which is 
characterised by the property that the action remains invariant under most small variations 
of a configuration. 
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The continuum limit is attained at /3 —>■ oo or 5 —)• 0. 

Our simulations were all performed with a cluster algorithm [TB] . Due to 
its non-local update steps, the Markov chain changes the topological sector 
frequently, which enables a precise direct measurement of the topological 
susceptibility. This result is then confronted with the value determined by 
the slab method, described in Section 2. 

We begin with an illustration of measured data for the probability product 
in eq. fl2.2p . as a function of the relative slab width x (here the “slab” is 
actually just an interval in Euclidean time). 

As a first example. Figured] shows data for (g^), as a function of x. These 
values were measured with the standard action at size L = 400 and /3 = 4, 
which implies a correlation length of ^ ~ 6.8. We see that the curves obtained 
in the sectors Q = 0, \Q\ = 1 and |Q| = 2 do follow the expected parabolas, 
which interpolate {q^){x = 0) = 0 and (g^)(x = 1) = as predicted in 
Section 2. 


Id 0(2), Standard Action 
L=400, |3=4 



Figure 1: The expectation value (g^)(a;), measured for the standard action 
at L = 400, /3 = 4, in the sectors |Q| = 0, 1, 2. We conhrm the predicted 
parabolic shape behaviour. 

Next we give results obtained with the constraint action, now at L = 100 
and 5 = 27r/3, where the correlation length is short, ^ ~ 1.1. Figure |2] shows 
again the numerical data for (g^) as a function of x, in the sectors (^ = 0, 
|(5| = 1 and IQ I = 2. For comparison, we also include the parabola for the 
function xtDx(l — x), which is predicted for (g'^)(x). Here we insert the di¬ 
rectly measured value of Xt- We see in all cases accurate agreement between 
this prediction and the numerical data. 
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Id 0(2), Constraint Action 

V=100, Sector Q=0, 5 =27i/3 



Id 0(2), Constraint Action 

V=100, Sector IQI=1, 8 =27t/3 



Id 0(2), Constraint Action 

V=100, Sector IQI=2, 8 =27:73 



Figure 2: The expectation values {q^){x) and (g'^)(a:), measured for the 
constraint action at L = 100, 5 = in the sectors |Q| = 0, 1, 2. For 

comparison, the black curve shows the predicted behaviour of (g'^)(a;), which 
agrees very well with the data. 

Let us proceed to a quantitative discussion. We start again with the stan- 
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dard action, where we consider (3 = 2 and /9 = 4. In a variety of sizes L we 
applied the slab method, and evaluated Xt by a £t of the data for to 

the predicted parabola. The results are given in Table [H Figure |3] shows the 
corresponding convergence of the scaling term Xt^, where ^ is the correlation 
length]^ They are consistently close to the analytical value at inhnite L (the 
corresponding formula is given in Ref. 0). That value is in all cases com¬ 
patible with the directly measured Xt, which is also included in Tabled! For 
increasing L the agreement with the slab method results improves further. 


L 

direct 

Q = 0 

IQI = i 

IQI = 2 

= 2 Xt,L=oo = 0.019364 ^ = 2.778866 

100 

0.019369(6) 

0.019280(7) 

0.019324(11) 

0.019480(29) 

200 

0.019372(8) 

0.019346(14) 

0.019314(10) 

0.019326(19) 

400 

0.019365(10) 

0.019348(16) 

0.019337(11) 

0.019333(19) 

13 = 4: Xt,L=oo = 0.007554 ^ = 6.814998 

150 

0.007554(3) 

0.007534(4) 

0.007550(6) 

0.007587(2) 

200 

0.007557(3) 

0.007541(3) 

0.007547(4) 

0.007581(8) 

250 

0.007549(3) 

0.007542(3) 

0.007550(4) 

0.007585(8) 

300 

0.007560(4) 

0.007545(3) 

0.007545(5) 

0.007554(5) 

350 

0.007554(5) 

0.007545(4) 

0.007557(2) 

0.007553(5) 

400 

0.007549(5) 

0.007552(4) 

0.007551(3) 

0.007564(5) 


Table 1: Explicit results for Xt by the slab method, for the standard action 
aX (3 = 2 and (3 = 4, as illustrated in Figure El 

Next we give analogous results for the Manton action. Figure El and Ta¬ 
ble [2| and for the constraint action. Figure [5] and Table El Qualitatively the 
same features are conhrmed. Quantitatively we see that the Manton action 
— which is classically perfect [9] — performs very well regarding the conver¬ 
gence towards the value at L = oo. For all actions, the convergence is best 
for IQ I < 1, while |Q| = 2 is affected by somewhat stronger hnite size effects. 


®Here and in Figures [2l and Owe insert the analytic values for which are also given 
in Tables m [21 and El They were calculated with the formulae given in Refs. [9] and [3] as 
functions of /3 and S, see also Tables mm and El This refers to L —>■ oo, but since L ^ ^ 
holds in all cases, this is not problematic. 
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1d 0(2), Standard action, p=2 
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Id 0(2), Standard action, p=4 
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Figure 3: The results for Xt, hi units of 1/^, based on the slab method, for 
the standard action at /3 = 2 (above) and /9 = 4 (below). For an increasing 
ratio L/^ they move towards the analytic result in the thermodynamic limit 
(dashed line). 
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We also verihed that (Q) and the kurtosis term (which is related to the 
Binder cumulant), 

C4 = 4(3{Qy-(Q‘‘)) , (3.3) 

are both compatible with 0 in all cases (within at most 2cr)I§ The former 


®For these parameters and huge statistics of 0(10®) measurements, the error in C 4 can 
be reduced to 0(10“®), and one observes significant deviations from zero [ 6 ]. However, 












Id 0(2), Manton action, p=2 
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Figure 4: Results for Xt, in units of 1/^, based on the slab method, for the 
Manton action at /3 = 2 (above) and /3 = 4 (below). Compared to the 
standard action, a signihcantly smaller ratio L/^ is sufficient for a very good 
approximation to the inhnite size result (dashed line). 
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follows from parity symmetry, and a Gaussian Q-distribution implies C 4 = 0. 
The vanishing of these two quantities is an assumption of the slab method. 
For the kurtosis this condition is not trivial; it will be considered in detail in 
the next section. 

Another source of systematic errors are artifacts in xt due to the hnite 
size and hnite lattice spacing. These artifacts are discussed in Appendix A. 

these tiny C 4 values have no influence on the interpretation of the results presented here. 
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L 

direct 

Q = 0 

IQI = i 

IQI = 2 

= 2 Xt,L=oo = 0.012663 ^ = 4.000321 

100 

0.012658(7) 

0.012666(5) 

0.012666(9) 

0.012631(22) 

200 

0.012661(4) 

0.012656(7) 

0.012661(6) 

0.012659(14) 

400 

0.012653(4) 

0.012672(7) 

0.012658(7) 

0.012666(11) 

13 = 4: Xt,L=oo = 0.006333 ^ = 8.000000 

150 

0.006330(4) 

0.006331(3) 

0.006332(4) 

0.006316(2) 

200 

0.006333(3) 

0.006330(2) 

0.006328(4) 

0.006339(9) 

250 

0.006335(3) 

0.006336(3) 

0.006333(3) 

0.006339(8) 

300 

0.006332(3) 

0.006333(4) 

0.006333(2) 

0.006316(9) 

350 

0.006334(4) 

0.006329(3) 

0.006335(4) 

0.006342(7) 

400 

0.006329(3) 

0.006330(3) 

0.006334(3) 

0.006323(5) 


Table 2: Explicit results for Xt by the slab method, for the Manton action at 
(5 = 2 and /3 = 4, as illustrated in Figure HI 


L 

direct 

Q = 0 

IQI = i 

IQI = 2 

5 = 27r/3 Xt,L=oo = 0.037037 ^ = 1.132367 

100 

0.037036(13) 

0.037281(16) 

0.037183(14) 

0.036990(35) 

200 

0.037008(15) 

0.037163(20) 

0.037121(18) 

0.037079(29) 

400 

0.037042(19) 

0.037091(33) 

0.037109(23) 

0.037056(26) 

5 = 1.5 Xt,L=oo = 0.018998 ^ = 2.451141 

100 

0.018987(8) 

0.019113(8) 

0.019058(10) 

0.018864(28) 

200 

0.019011(9) 

0.019050(9) 

0.019049(13) 

0.018983(17) 

400 

0.018992(7) 

0.019038(13) 

0.019032(8) 

0.019019(12) 

5 = 1 Xt,L=oo = 0.008443 ^ = 5.793617 

100 

0.008443(5) 

0.008491(3) 

0.008416(10) 

0.008208(36) 

200 

0.008445(4) 

0.008475(5) 

0.008452(7) 

0.008397(16) 

400 

0.008439(5) 

0.008460(5) 

0.008455(6) 

0.008445(15) 


Table 3: Explicit results for Xt by the slab method, for the constraint action 
at (5 = 27r/3, 1.5 and 1, as illustrated in Figure [5l 
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1d 0(2), Constraint action, 5=27t/3 
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Figure 5: Results for Xt, in units of 1/.^, based on the slab method, for the 
constraint action at 5 = 27r/3, 1.5 and 1, as a function of L/^. 
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4 Results for the 2d Heisenberg model 

We proceed to the 2d 0(3) model, where we consider square lattices with 
classical spins G on the sites x. In order to dehne the topological charge, 
we cut each plaquette into two triangles (with an alternating orientation 
between adjacent plaquettes). For a triangle with sites x, y, z we identify 
the (minimal) solid angle spanned by Cx, Cy, Cz, including a sign factor for 
its orientation (a fully explicit description is given in Ref. i)- 

The sum of the two oriented solid angles within a plaquette, divided by 
47 r, dehnes the topological charge density Qx- The total charge of a spin con- 
hguration [e] amounts to Q[e] = ^ counts how many times the 

sum of these solid angles covers the sphere 8“^ with a dehnite orientation. 

In analogy to Section 3, we consider three lattice actions, including an 
obvious generalisation of the Manton action, which was simulated with a new 
variant of the Wolff cluster algorithm, 

'y ^ (1 ' (^x+ji) 

x,fi=l,2 

/3 ^ (l - [arccos(ea, • Cx+f,)]'^'^ 

f 0 Cx- ex+ji > cos 5 Vx, y = l,2 
[ +00 otherwise. ^ ' 

Kegaiding the application of the slab method with these lattice actions. 
Figure [6] shows results obtained in the lattice volumes V = 48^ and 64^. 

We also consider further volumes, including rectangular shapes. In fact, 
it is not obvious if the slab method results should be compared to Xt in the 
entire (periodic) volume V, or to Xt,v /2 in a (non-periodic) slab of size V/2. 
Hence we have measured the latter as well, for comparison. The detailed 
results, alon^ with the corresponding correlation length, are given in Tables 
m |5] and |6][] However, for volumes, which are large enough for the slab 
method to work quite well, it turns out that Xt and Xtyji are too close to 
each other to be distinguished in light of the slab method results. 

The assumption of a Gaussian distribution of the topological charges is 
essential for the viability of this method. Figure [7] shows that the kurto- 

^In this model, the apparent scaling quantity Xt diverges logarithmically in the 
continuum limit (see e.g. Ref. [3]), so here we refer directly to Xt- 


‘S'standard [t- ] 

‘S'Manton [c ] 

‘^'constraint \c ] 
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Figure 6: Results for the topological susceptibility Xt in the 2d 0(3) model, 
in the volumes V = 48^ and 64^, for three lattice actions. Each set of 
results consists of four data points and shows (from left to right) the directly 
measured value, and the results by the slab method in the sectors \Q\ = 
0 , 1 , 2 . 

sis fl3.3p — as a measure for the deviation from a Gaussian — decreases 
rapidly as we approach the continuum limit; for the Manton action we see 
an exponential decrease of |c 4 | as increases. Thus a non-zero value of |c 4 | 
appears as a lattice artifact, and not as a hnite size effect; this observation 
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/S 

R 

^direct 

Xt 

^direct 

Xt,\//2 

Q=o 

Xt 

Xt 

— 

1 

48^ 

0.01715(2) 

0.01718(2) 

0.01714(2) 

0.01723(2) 

0.01705(2) 


642 

0.01716(3) 

0.01717(2) 

0.01718(2) 

0.01717(2) 

0.01726(2) 


642 

0.002319(3) 

0.002354(2) 

0.002305(5) 

0.002322(5) 

0.002361(5) 


802 

0.002332(5) 

0.002356(3) 

0.002327(4) 

0.002327(4) 

0.002358(4) 

1 ^ 

64 X 128 

0.002332(4) 

0.002344(3) 

0.002304(3) 

0.002321(2) 

0.002324(3) 


962 

0.002341(3) 

0.002359(3) 

0.002333(3) 

0.002337(3) 

0.002352(3) 


1282 

0.002327(3) 

0.002347(3) 

0.002390(10) 

0.002381(9) 

0.002371(9) 


128 X 256 

0.002334(4) 

0.002342(3) 

0.002353(1) 

0.002351(1) 

0.002366(1) 


Table 4: Results for the topological susceptibility Xt in fhe 2d 0(3) model 
with the standard action. We consider a variety of volumes at /9 = 1 and 
at /5 = 1.5, with a correlation length of ^ = 1.3 and 9.4, respectively. The 
results for V = 48^ and 64^ are illustrated in Figure |6l 


/S 

R 

^direct 

Xt 

^direct 

Xt,\//2 


Wl=i 

Xt 

|y|=2 

Xt 

0.5 

48^ 

64^ 

0.02329(4) 

0.02335(2) 

0.02333(3) 

0.02335(3) 

0.02330(3) 

0.02330(2) 

0.02339(2) 

0.02343(2) 

0.02330(2) 

0.02330(2) 

0.75 

48^ 

64^ 

0.00954(2) 

0.00953(2) 

0.00956(1) 

0.00957(1) 

0.00953(2) 

0.00962(2) 

0.00952(2) 

0.00955(1) 

0.00952(2) 

0.00957(2) 


Table 5: Results for the topological susceptibility Xt in the 2d 0(3) model 
with the Manton action. We consider two volumes at /3 = 0.5 and at (3 = 0.75, 
with a correlation length of ^ = 0.9 and 1.7, respectively. These results are 
illustrated in Figure |6l 


agrees qualitatively with a discussion in Ref. [B]. This artifact is much less 
suppressed for the standard action and the constraint action. Regarding the 
application of the slab method with these lattice actions. Figure [6] shows re¬ 
sults obtained in the lattice volume V = 48^ and 64^. The systematic errors 
of the slab method appear as tiny deviations of Xt extracted from |Q| =0, 1 
and 2 (three rightmost points in each set of four) from the directly measured 
value (leftmost point). 
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5 

R 

^direct 

Xt 

^direct 

Xt,\//2 


Xt 

— 

Ptt 

48^ 

0.02157(5) 

0.02155(2) 

0.02145(2) 

0.02162(2) 

0.02149(2) 

!{ j Kj 

642 

0.02144(4) 

0.02151(3) 

0.02144(2) 

0.02149(2) 

0.02156(2) 

9 

482 

0.01545(2) 

0.01547(2) 

0.01549(2) 

0.01546(1) 

0.01543(2) 


642 

0.01549(3) 

0.01547(2) 

0.01538(2) 

0.01541(1) 

0.01543(1) 


162 

0.002445(6) 

0.002663(4) 

0.002168(65) 

0.003806(14) 

0.005211(27) 


322 

0.002795(4) 

0.002873(3) 

0.002816(17) 

0.002888(16) 

0.003146(15) 

n TT 

642 

0.002797(4) 

0.002835(4) 

0.002859(6) 

0.002853(7) 

0.002850(7) 


64 X 128 

0.002795(4) 

0.002816(3) 

0.002826(3) 

0.002812(3) 

0.002810(3) 


962 

0.002792(5) 

0.002818(3) 

0.002831(4) 

0.002835(4) 

0.002843(4) 


1282 

0.002783(6) 

0.002805(3) 

0.002837(3) 

0.002828(3) 

0.002819(3) 


Table 6: Results for the topological susceptibility Xt in the 2d 0(3) model 
with the constraint action. We consider a variety of volumes at 5 = 27r/3, 2 
and 0.557r, with a correlation length of .^ = 0.85, 1.05 and 3.5, respectively. 
The results for V = 48^ and 64^ are illustrated in Figure |6l 
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Figure 7: The absolute value of the kurtosis, given in eq. (13.3p . in the 2d 
0(3) model. For each of the three actions, and at a hxed correlation length 
we display the results mV = 48^ (on the left), and V = 64^ (on the right). 
We see that the volume hardly affects the value of |c 4 |, which decreases for 
increasing correlation length This convergence to zero is fastest for the 
Manton action, where we recognise an exponential decrease. 
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5 Conclusions 


We have tested an unconventional method for the numerical measurement 
of the topological susceptibility Xt, based on the division of the volume into 
slabs. This method is applicable even if conhgurations in only one topological 
sector are available. 

Our study shows that — under suitable conditions — this method works 
very well. Its statistical accuracy is comparable to the precision in the ab¬ 
sence of “topological slowing down”, i.e. the slab method is not affected by 
the freezing of the topological sectors. 

In particular, we obtained in the 2d 0(3) model results for Xt-, which are 
correct on the percent level, and in the Id 0(2) far beyond. This precision 
was revealed by large statistics of O(IO^) measurements, which would not be 
accessible in higher dimensional models like QCD; in this sense, the accuracy 
of the slab results are fully satisfactory. 

The slab method is most successful at small topological charges, IQI < 
1, whereas its application in the sector |(5| = 2 is more tedious, since the 
hnite size effects are worse. However, even at |Q| < 1 the hnite size effects 
are highly persistent. In fact, other methods and formulae show that these 
effects are only suppressed by a power series m 1/V for topologically hxed 
measurements OE]. 

As an illustrative example of possible systematic effects, we refer to the 
results at 5 = O.SStt in Table [6j even the smallest volume, V = 16^, cor¬ 
responds to the ratio L/^ ~ 4.6, where one would expect usual (i.e. expo¬ 
nentially suppressed) hnite size effects to be mostly eliminated. The results 
for the slab method, however, are strongly distorted — they improve signih- 
cantly at L = 32 (and for |(5| = 2 it even takes L = 64). 

This illustrates the general feature: at hxed (3 or 6, i.e. at approximately 
constant the slab results converge for increasing volume. In order to assure 
that they do converge to the (vicinity of the) correct value, we further have to 
require the topological charges Q to be (approximatively) Gauss-distributed 
(that also ensures a good ht to formula 12.2p . Therefore we also studied 
the kurtosis, which may deviate from zero (indicating a deviation from a 
Gaussian) as a lattice artifact. It is suppressed when ^ increases (in lattice 
units), with a rate that depends on the lattice action; for the Manton action 
we observed a particularly fast, exponential suppression. 

Hence ^ should be sufficiently large to control that requirement, and the 
volume should be large enough to obtain L/^ = 0(10); then the slab method 


16 


can be expected to provide the correct Xt value on the percent level. 
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A Artifacts in the topological susceptibility 

An essential point for the quality of the results obtained by the slab method 
is the precision of Xt as measured in the sub-volumes. In this appendix 
we discuss the artifacts that occur: Sub-Appendix A.l deals with the hnite 
size effects, which depend on the ratio and which even occur in the 
continuum formulation. Sub-Appendix A.2 considers the lattice artifacts, as 
a function oia/L (where a is the lattice spacing, and L is now the dimensional 
size). Both considerations are performed in the analytically tractable case of 
the quantum rotor; Sub-Appendix A.2 captures the standard and the Manton 
action. For simplicity we deal with periodic boundary conditions throughout, 
although they are partly open in the sub-volumes, so we assume implicitly 
that L is large enough for this distinction to be negligible. 

A.l Finite size effects 

While the slab method is plagued by particularly persistent hnite size effects, 
the true Xt value tends to converge exponentially as the volume is enlarged. 
Here we discuss this convergence for the case of the Id 0(2) model. 

In the continuum formulation, with the action S'):/?] = f /g 0^ dt, the 
topological susceptibility takes the form |9] 



a := 27r/5/L . (A.l) 


It can be written in terms of a Jacobi 0-function, 




(A.2) 
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By applying the Jacobi identity, 


e{a) = 


‘I'Jiyfa 


0(l/a) , 


(A.3) 


we obtain 


0 («)|c 




(^1 + 2e-"/“ + ... ) 


Inserting approximation flA.jp into formula (IA.2p . we arrive at 

1 L 


Xt ^ 


_ g-L/(2/3) ^ g-L//3 


(27r)2/5 2'k‘^I3‘^ 

With the correlation length ^ = 2/3 [9], this corresponds to the ratio 


Xt 


4L 


^ 1 + f _ 

Xt,L=oo S 


(A.4) 


(A.5) 


(A.6) 


Thus we see explicitly the exponentially suppressed hnite size corrections. 
Figure [8] illustrates this formula, and compares it to the exact result (nu¬ 
merical summation of the series in eq. flA.ip i. which agrees with simulation 
data for the Manton action at /3 = 4 (where lattice artifacts are practically 
erased). 


A.2 Lattice artifacts 


We begin with general properties of a held theory in a periodic Euclidean 
volume V, where the conhgurations of the held 0 are divided into topological 
sectors. In presence of a vacuum angle 0 (and with h = 1), the partition 
function and the topological susceptibility can be written a^ 


Z{e) = y ^ 


1 

r Z" /Z'\2i 


. z ^\z). 


(A.7) 


In the Hamilton formulation, the partition function reads 


® Occasionally we omit the argument 9 in the partition function and (below) in the 
energy eigenvalues. 
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Finite size effects of Xt in the 1d 0(2) model 



Figure 8: The ratio XtlXt,L=oo as a function of ^/L. We show the leading 
and next-to-leading hnite size correction in approximation flA.6j) . the exact 
result and simulations results for the Manton action at = 4. 


where (3 is the inverse temperature, and we assume a discrete energy spec¬ 
trum. Inserting eq. flA.7p leads to 

» = IvEK - + ■ (A-8) 

n n 

Now we focus on d = 1, and replace V and (3 both by L (the periodicity 
range in Euclidean time). More specihcally, we consider again the quantum 
rotor, where L also represents the moment of inertia. We deal with a lattice 
of spacing a, and periodicity over N sites, L = Na. A consideration of the 
transfer matrix yields [9] 



where n G Z, and 


fi^) = 


1 — cos (p 


standard action 
Manton action. 


In either case, f{(p) is even, such that E'^{0) = 0, in accordance with ((5)e=o = 
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0, and some algebra leads to 


I ^2l rAf-2 


Xt,0=O — 


iE„[-& + (iv + i)Km 


TN 

^ Z_-/m '^m 

/ TT 

cos(n(^) 


/^(O) = 


dip e 




dip e /(‘^)V“(^2 _ (A.9) 


To proceed, we now have to specify f{p). We start with the Manton 
action, and substitute p = p\fN, 


4(0) = 


^/N J-nPW 

-1 /‘tt'n/TV 


d(p e cos{n(f/\^) = — 


/;(0) = ^y_^^#^■"^'Vsi^(n^/^4v) = ^J« 

1 rn^N -I 


Manton _ 

At,6=0 — 


{2TifLN 




^(0)'|V 
m ) 


(A. 10) 


If we write cos(np/\fN) and sm.(np/\fN) in exponential form, and complete 
the squares in the exponents, 

-]^p^+ mpl\fN =-]^{p-m/\fNf - ^ , (A.ll) 

we see several types of lattice artifacts, such as the extra term —an? j (2L) in 
the exponent, and the incomplete Gauss integrals; they are exponentially sup¬ 
pressed in a/L. The substitution of the integration variable p ^ p — in/y/N 
corresponds to a shift of the integration contour; for the incomplete Gauss 
integral this is another artifact, which is exponentially suppressed. In sum¬ 
mary, there are no power lattice artifacts for the Manton action. 
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For the standard action the corresponding formulae read 


standard _ 

Xt,e=o — 


with 


1 + {N+ 

(27r)2a 

= f d(p e(“"^-i)^cos(n(p) 

J —TT 

/ TT 

dip e(“^‘"-^)^(psin(n(p) 

■TT 

j(2) = /" d(pe(“"^-i)Vcos(n(p) . 


(l)^2l/ r(0)\iV-2 


(A. 12) 


Regarding the search for power lattice artifacts, we start from the Manton 
action and add the two leading modihcations of the standard action, f{ip) — 
(p^/2 —(p"^/4! + (p®/6!. This modihes the expressions for /„(0), /^(O) and /"(O) 
in eq. flA.lOD by a factor 

+ 2 ( 40 ^ 5 *'* - 6 !^'"') 

under the integrals, and after some calculus we obtain 

xtsr" ^ xttT(l + + (A,13) 


Hence for the standard action Xt, 9 =o is affected by corrections linear in a/L. 

This may appear surprising for a bosonic theory, but one has to keep in 
mind that Xt, 9 =o is not a scaling quantity. The actual scaling artifacts refer 
to the product XtC fh® thermodynamic limit they are expressed in powers 
of a/f |3], 


xti = < 


V 


1 

27r^ 


^ + ||j ... j standard action 

y/Tifa{l — + ... j Manton action 

^ . j constraint action. 


For the standard action, the linear artifacts cancel in this scaling quantity, 
while the Manton action is classically perfect; its scaling artifacts are expo¬ 
nentially suppressed. As an exotic feature, the constraint action does have 
linear artifacts, due to the absence of a derivative term in the action. 


We summarise that the hnite size effects are exponentially suppressed in 
L/^, whereas the type of lattice artifacts depend on the lattice action. 


21 










References 


[1] E. Witten, Nucl. Phys. B 156 (1979) 269. G. Veneziano, Nucl. Phys. B 
159 (1979) 213. 

[2] M. Liischer, Nucl. Phys. B 549 (1999) 295 [hep-lat/9811032]. 

[3] W. Bietenholz, U. Gerber, M. Pepe and U.-J. Wiese, JHEP 12 (2010) 
020 [arXiv: 1009.2146 [hep-lat]]. 

[4] O. Akerlund and P. de Forcrand JHEP 1506 (2015) 183 

[arXiv:1505.02666 [hep-lat]]. 

[5] R. Brower, S. Ghandrasekharan, J.W. Negele and U.-J. Wiese, Phys. 
Lett. B 560 (2003) 64 [hep-lat/0302005]. 

S. Aoki, H. Fukaya, S. Hashimoto and T. Onogi, Phys. Rev. D 76 (2007) 
054508 [arXiv:0707.0396 [hep-lat]]. 

S. Aoki et al. (JLQGD and TWQGD Gollaborations), Phys. Lett. B 665 
(2008) 294 [arXiv:0710.1130 [hep-lat]]. 

W. Bietenholz, 1. Hip, S. Shcheredin and J. Volkholz, Eur. Phys. J. C 
72 (2012) 1938 [arXiv;1109.2649 [hep-lat]]. 

A. Dromard and M. Wagner, Phys. Rev. D 90 (2014) 074505 
[arXiv: 1404.0247 [hep-lat]]. 

H. Fukaya, S. Aoki, G. Gossu, S. Hashimoto, T. Kaneko and J. Noaki 
(JLQGD Gollaboration), arXiv:1411.1473 [hep-lat]. 

[6] 1. Bautista, W. Bietenholz, A. Dromard, U. Gerber, G.P. Hofmann, H. 
Mejia-Diaz and M. Wagner, arXiv:1503.06853 [hep-lat]. 

[7] P. de Forcrand, M. Garcia Perez, J.E. Hetrick, E. Laermann, J.F. La- 
gae and 1.0. Stamatescu, Nucl. Phys. Proc. Suppl. 73 (1999) 578 [hep- 
lat/9810033]. 

[8] R.G. Brower et al. (LSD Gollaboration), Phys. Rev. D 90 (2014) 014503 
[arXiv:1403.2761 [hep-lat]]. 

[9] W. Bietenholz, R. Brower, S. Ghandrasekharan and U.-J. Wiese, Phys. 
Lett. B 407 (1997) 283 [hep-lat/9704015]. 

[10] S. Diirr, Z. Fodor, G. Hoelbling and T. Kurth, JHEP 0704 (2007) 055 
[hep-lat/0612021]. 


22 



[11] R. Sinclair, Phys. Rev. D 45 (1992) 2098. 

[12] B. Berg and M. Liischer, Nucl. Phys. B 190 (1981) 412. 

[13] N. Manton, Phys. Lett. B 96 (1980) 328. 

[14] W. Bietenholz, M. Bogli, F. Niedermayer, M. Pepe, F.G. Rejon-Barrera 
and U.J. Wiese, JP[EP 1303 (2013) 141 [arXiv: 1212.0579 [hep-lat]]. 

W. Bietenholz, U. Gerber and F.G. Rejon-Barrera, J. Stat. Mech. (2013) 
P12009 [arXiv: 1307.0485 [hep-lat]]. 

[15] U. Wolff, Phys. Rev. Lett. 62 (1989) 361. 


23 



